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Abstract. The global solutions in critical spaces to the multi-dimensional compressible 
viscoelastic flows are considered. The global existence of the Cauchy problem with 
initial data close to an equilibrium state is established in Besov spaces. Using uniform 
estimates for a hyperbolic-parabolic linear system with convection terms, we prove the 
global existence in the Besov space which is invariant with respect to the scaling of the 
associated equations. Several important estimates are achieved, including a smoothing 
effect on the velocity, and the L 1 — decay of the density and deformation gradient. 



1. Introduction 

We consider the following equations of multi-dimensional compressible viscoelastic flows 

Dana ED: 

Pt + div(pu) = 0, (1.1a) 

(pu)t + div (pu ® u) - pAu- (A + p)Vdivu + VP(p) = adiv(pFF T ), (1.1b) 

Ft + u-VF = VuF, (1.1c) 

where p stands for the density, u G M. N (N = 2, 3) the velocity, and F 6 M NxN (the set of 
N x N matrices) the deformation gradient. The viscosity coefficients p, A are two constants 
satisfying p > 0, 2p + iVA > 0, which ensures that the operator — /iAu — (A + p) Vdivu is a 
strongly elliptic operator. The pressure term P(p) is an increasing and convex function of 
p for p > 0. The symbol (S> denotes the Kronecker tensor product, F T means the transpose 
matrix of F, and the notation u • VF is understood to be (u • V)F. For system (|1.1|) . 
the corresponding elastic energy is chosen to be the special form of the Hookean linear 
elasticity: 

W(F) = ^\F\ 2 , a>0, 

which, however, does not reduce the essential difficulties for analysis. The methods and 
results of this paper can be applied to more general cases. 

In this paper, we consider the Cauchy problem of system (II. ip subject to the initial 
condition: 

(p,u,F)| t=0 = OSo(x),u (x),F (x)), xeR N , (1.2) 

and we are interested in the global existence and uniqueness of strong solution to the 
initial- value problem (jl.ip - (jl.2p near its equilibrium state in the multi-dimensional space 
M. N . Here the equilibrium state of the system (jl.ip is defined as: p is a positive constant 
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(for simplicity, p = 1), u = 0, and F = I (the identity matrix in M 3x3 ). We introduce a 
new unknown variable E by setting 

F = I + E. 

Then, (jl.ip becomes 

pt + div(pu) = 0, (1.3a) 

(pu)t + div (pu <g) u) - fiAu - (p + A)Vdivu + VP(p) = adiv(p(7 + £;)(/ + £) T ), 

(1.3b) 

Et + u- VE = VuE + Vu, (1.3c) 
with the initial data 

(p,u,E)\ t=0 = (p (x),u (x),E (x)), xeR N . (1.4) 

There have been some results about the local existence of strong solutions to the com- 
pressible viscoelastic flows, see |10| [T5] and the references therein. The global existence to 
(II. ip is a difficult problem due to the appearance of the deformation gradient. The chal- 
lenge is to identify an appropriate functional space where the Cauchy problem (|l.ip - (|1.2p 
is well-posed globally in time. In this paper, to construct a global solution, we are going to 
use the scaling for the compressible viscoelastic flow to guess which space may be critical. 
We observe that system (jl.ip is invariant under the transformation 

(p (x),u (x),F (x)) (p (lx),lu (lx),F (lx)), 
(p(t,x),u(t,x),F(t,x)) -)• (p(l 2 t,lx),lu(l 2 t,lx),F(l 2 t,lx)), 
up to changes of the pressure law P into l 2 P, and a into l 2 a. This suggests the follow- 
ing definition: A functional space 21 C S'(R N ) x (S'(R N )) N x (S'(R N )) NxN is called a 
critical space if the associated norm is invariant under the transformation (p, u, F) — >• 
(p(l-),lu(l-),F(l-)) (up to a constant independent of I), where S' is the space of tem- 
pered distributions, i.e., the dual of the Schwartz space S. According to this definition, 
flf x (fif- 1 )" x ST (see Section 2 for the definition of B s := B^ X (R N )) is a criti- 
cal space. The motivations to use the homogeneous Besov space B s with the derivative 
index y include two points: first, is an algebra embedded in L°°, which allows us 
to control the density and the deformation gradient from below and from above without 
requiring more regularity on derivatives of p and F; second, the product is continuous from 
BT ~ a x £T to Bi ~ Q for < a < N. 

For the global existence, the hardest part of the argument is to deal with the linear 
terms Vjo, divE and Vu, especially the first two terms. It turns out that finding some 
dissipation for divE is a crucial step. This step for the incompressible case has been 
fulfilled successfully in [H]. For the compressible system (jl.ip . we will reformulate the 
system, and use the divergence-free property of compressible viscoelastic flows for the 
"compressible" part of the velocity, while the property on curl is used to deal with the 
"incompressible" part of the velocity. Meanwhile, we decompose the deformation gradient 
into two parts: the symmetric part and the antisymmetric part. With this technique 
and decomposition, we will be able to obtain successfully the dissipation estimates on the 
density and the deformation gradient for an auxiliary system with convection terms. These 
estimates are crucial for the global existence. We remark that for the global existence of 
solutions to (II. lh near equilibrium, the intrinsic properties of the divergence and the curl 
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are important and necessary. During the final stage of this paper, we noticed that some 
similar results are also obtained independently in [20], where the intrinsic properties of 
the divergence and curl of viscoelastic flows (see Appendix) to system (II. ip are used to 
control the dissipation of the deformation gradient F. 

For the incompressible viscoelastic flows and related models, there are many papers 
in literature on classical solutions (cf. [7J [121 031 E] an d the references therein). On 
the other hand, the global existence of weak solutions to the incompressible viscoelastic 
flows with large initial data is still an outstanding open question, although there are 
some progress in that direction ([TBI f!8l 119]). For the well-posedness of global solutions 
to the compressible Navier-Stokes equations, see [31 0] (for barotropic cases), and [5] (for 
barotropic cases with heat conduction). For the inviscid elastodynamics, see [22] and their 
references on the global existence of classical solutions. 

The rest of this paper is organized as follows. In Section 2, we review the definitions 
of Besov spaces and show some good property of the Besov spaces. In Section 3, we 
reformulate the system (jl.ip and state the main theorem. Section 4 is devoted to a priori 
estimates for an auxiliary linear system with convection terms. In Section 5, we give 
the proof of our main result, while in the appendix (Section 6), we prove two intrinsic 
properties of compressible viscoelastic flows. 



2. Basic Properties of Besov Spaces 

Throughout this paper, we use C for a generic constant, and denote A < CB by A < B. 
The notation A ~ B means that A < B and B < A. Also we use (a q ) q ^z to denote a 
sequence such that X^ez Q <2 — ^' (/Iff) denotes the inner product of two functions /, g 
in L 2 (R Ar ). The standard summation notation over the repeated index is adopted in this 
paper. 

The definition of homogeneous Besov spaces is built on an homogeneous Littlewood- 
Paley decomposition. First, we introduce a function ip £ C°°(M. N ), supported in the shell 

c = U€i^<iei<f}, 

such that 

;>>(2- 9 o = i, ^ e t^o. 

Denoting by h := J 7 ~ 1 ip the inverse Fourier transform of ip, we define the dyadic blocks 
as follows: 

A q f = rJ>(2~*D)f = 2* N [ h(2Py)f(x - y)dy, 

and 

sj= ^2 a p/' 

p<q-l 

where D is the first order differential operator. The formal decomposition 



/ = £ A «/ 



(2.1) 
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is called the homogeneous Littlewood-Paley decomposition. But unfortunately, the above 
identity is not always true in S'(R N ) as pointed out in @]. Nevertheless, (|2.ip is true 
modulo polynimials (see [H El [6] ) . 

For s£l and / G S'(R N ), we denote 

B s :=Y,^\\\f\\ L ,. 



Notice that || • \\b* is only a semi-norm on {/ G S'{U. N ) : ||/||b s < because ||/||b s 
vanishes if and only if / is a polynomial. This leads us to introduce the following definition 
for homogeneous Besov spaces: 

Definition 2.1. Let s£l and m = — [y + 1 — s] . If m < 0, we set 



B s 



f G S'(R N ) : \\f\\ B s < oo and / = in S'( 



If m > 0, we denote by V m the set of polynomials with N variables of degree < m and 
define 



f G S'(R N )/V m : ||/||b. < oo and / = ]T A,/ in S'(R N )/V n 



Functions in -B s have many good properties (see Proposition 2.5 in [4]): 

Proposition 2.1. T/ie following properties hold: 

• Density: the set is dense in B s if \s\ < y; 

• Derivation: ||/||b« ~ ||V/||b=-i/ 

• Fractional derivation: let T = \J — A and <r G R; i/ien i/ie operator r°" is an 
isomorphism from B s to B s ~ a ; 

• Algebraic properties: for s > 0, B s D L°° is an algebra; 

• Interpolation: (B Sl , B S2 ) e>1 = B fisi+ ( 1_ ^ Sa . 

For the composition in we refer to [2] for the proof of the following estimates: 

Lemma 2.1. Giuen s > and f G L°° n 5 s . 

• Let^ e W^ C +2 (R N ) such that <5(0) = 0. T/ien *(/) G 5 s . Moreover, there exists 
a function C of one variable depending only on s, N and ^f, and such that 

||*(/)||B.<C(||/|| £ »)||/|| fl .. 

• Let <I> G Wj° ] c +2 (R N ) such that $'(0) = 0. Suppose that f and g belong to B^ and 
that (/ - g) G B s for some s G (— y, y]- Then <£(/) - $(5) belongs to B s and 
there exists a function of two variables C depending only on s, N and <!>, and such 
that 

\Mf) - <&( 5 )|| B , < C(\\f\\ L «,, l^llico) (\\f\\ B M + ) 11/ - g\\ B s. 

But, different from the nonhomogeneous Besov space, the homogeneous Besov spaces 
fail to have nice inclusion properties. For example, owing to the low frequencies, the 
inclusion B s B r does not hold for s > r. Still, the functions of B s are locally more 
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regular than those of B r : for any ip £ Cq° and / S B s , the function iff is in B r . This 
motivates the definition of hybrid Besov spaces where the growth conditions satisfied by 
the dyadic blocks are not the same for low and high frequencies. Let us recall that using 
hybrid Besov spaces has been crucial for proving global well-posedness for compressible 
gases in critical spaces (see [U [5] ) . The definition of the hybird Besov space is given as 
follows (see Definition 2.8 in [3] or [5]). 

Definition 2.2. Let s, t 6 R. We set 

ii/ii^ = E 29S ii a ^ii^+E 29< ii a ^ii^- 

<?<o g>0 
Denoting m = — [y + l — s] , we define 

B s * = {/ € S'(R N ) : ||/||a., t < oo} if m < 0, 
B^ = {/ € S'(R N )/V m : II/Hb.,* < oo} if m > 0, 

Remark 2.1. Some remarks about the hybrid Besov spaces are in order: 

• B s ' s = B s ; 

• Ifs < t, thenB s '< = B s nB*. Otherwise, B s ' 1 = B s +B l . In particular, B s >f ^ L°° 

clS 5 ^ j 

• The space B 0,s coincides with the usual nonhomogeneous Besov space 



/ G S'(R N ) : \\ X (D)f\\ L 2 + £VlA,/|| i2 < oo i , where = 1 - ^0(2^); 

9>0 J <2>0 

• If si < s 2 and t x > t 2 , then B 81 '* 1 ■-)• B S2 '* 2 . 

For products of functions in hybrid Besov spaces, we have (see Proposition 2.10 in [4]): 

Proposition 2.2. Given s±, S2,ti,t 2 £ H^- 

• For a// s\,S2 > 0, 

ll/^bn^ ^ II/IU°°IMIb>i,s 2 + IMU°°ll/llB«i,»a- 

• For a// s±,S2 < y such that minjsi + 1%, S2 + t 2 } > 0, 

\\fg\\g sl+S2 -f ,i 1+ t 2 -f ^ 11/11 ^-1,-2 ll^ll s*i .*2 ■ 

In order to state our existence result, we introduce some functional spaces and explain 
the notations. Let T > 0, r G [0, oo] and X be a Banach space. We denote by A4(0, T; X) 
the set of measurable functions on (0, T) valued in X. For / € M(0,T; X), we define 

rT 



mx) = (J WfirWxdrY ifr<oo, 
=supess r6(0j T)||/(T)||x- 



Denote 

L r (0,T;X) = {/ G X(0,T;X) : ||/|| L r (x) < oo}. 

If T = oo, we denote by L r (M. + ;X) and ||/||i,rpn the corresponding spaces and norms. 
Also denote by C([Q,T],X) (or C(M. + ,X)) the set of continuous X-valued functions on 
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[0, T] (resp. R + ). We shall further denote by Cfe(]R + ;X) the set of bounded continuous 
X- valued functions. 

For a e (0,1), C a ([0,T];X) (or C a (R + \X)) stands for the space of the Holder con- 
tinuous functions in time with order a, that is, for every t,s in [0,T] (resp. we 
have 

\\f(t)-f(s)\\x<\t-s\ a . 

In this paper, the following estimates for the convection terms arising in the localized 
system will be used several times (cf. Lemma 5.1 in [5] or Lemma 6.2 in [3]). 



Lemma 2.2. Let G be an homogeneous smooth function of degree m. Suppose — ^ < 

N 
2 

\(G(D)A q (u.Vf)\G(D)A q f)\ 



Si,U < 1 + for i = 1,2. Then the following inequalities hold: 



<Ca ? 2-^ sl,S2 W— )|| u || Dl+4 ||/||^ liS2 \\G(D)A q f)\\ L2 , 



(2.2) 



and 



where 



\(G(D)A q (u ■ Vf)\A q g) + (A,(u • Vg)\G(D)A q f)\ 

< Ca q \\u\\ Bl+ u (2^' i2 W- m ) \\G(D)A q f\\ L2 \\g\\ §H , t2 (2 3) 



s, if q< 0, 
t, if q>l. 



For the nonlinear term VuE, we have the following estimates. 
Lemma 2.3. // -f < st < 1 + f for i = 1, 2, then 

\\VuE\\ Sai ,s 2 < C\\u\\ b1+ n\\E\\ §s1 , S2 . (2.4) 

To proof the above lemma, we need to recall the paradifferential calculus which enables 
us to define a generalized product between distributions. The paraproduct between / and 
g is defined by 

T f9 = ^2S q - 1 fA q g. 

We also have the following formal decomposition (modulo a polynomial): 

fg = T f g + T g f + R(f,g), 

with 



where A q = A q _i + A fy + A g+1 . 

Proof of Lemma \2. 3[ Denoting T'j-g = Tfg + R(g, /), we get the following decomposition 

A q (VuE) = A q T' E Vu + J q , (2.5) 
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with 

J q = Yl ([A„^„ 1 (Vu)]A g ^ + (5 g ,_ 1 -5 g _ 1 )(Vu)A g A (? ^ + 5 (? _ 1 (Vu)A^). 

|g'-«|<3 

Applying Proposition 5.2 in [5] or Proposition 6.1 in [3], we see that the first term on the 
right-hand side of ()2.5[) satisfies (|2.4p provided — y <«i < y + ^ ^ or * = 1> 2- Next we 
estimate J q term by term. First, for the commutator [A q , S q i-i(Vu)]A q /E, we have 

[A q ,S ql ^(Vu)]A ql E(x) 

= 2"*/ [ h{y)(yS q ^ l (VVu)(x-2- (1 Ty))A q ,E{x-2-iy)dTdy. 

JR N JO 

The above identity, together with Young's inequality for convolution operator, yields 
||[A ? ,S 3 >-i(Vu)]A 3 ,£;|| L2 < CUVullLooHAg/^llia, 

since 

||^/_iWu||loo < 2< z ||Vu|| i ° < 2 q \\u\\ n +1 

according to Bernstein's Lemma (cf. [HE]). Hence, we easily obtain the following inequal- 
ity: 

||J 9 || i2 <C« g 2-^ sl ' S2 W||u|| sl+f \\E\\ §S1 , S2 . 
The proof is complete. □ 



3. Reformulation and Main Results 

In this section, we state our global existence result. We first reformulate system f) 1 . 1 j) . 
Assume that the pressure P(p) is an increasing convex function with P'(l) > 0, and denote 
Xo = For p > 0, system (jl.ip can be rewritten as 

pt + u • Vp + pdivu = 0, (3.1a) 

d t Ui + u • Vtii - i (pAui - (A + p)d Xi dwu) + —^-d Xi p = aF jk d x .F ik , (3.1b) 

P P 
F t + u-VF = VuF, (3.1c) 

where we used the condition div(pF T ) = (see Lemma l6.ip for all t > 0, which ensures 
that the i-th component of the vector div(/oFF T ) is 

dxjipFikFjk) = pFjkd Xj Fik + ?ikd Xj (pFjk) 
= pPjkd Xj Fik- 

Define 

p(t,x) = p(xlt,Xox) - 1, u(t,x) = Xou(xo*>XoaO) = F (Xo t > Xo^) - J, 

then 

pt + u • Vp + divu = — pdivu, (3.2a) 
a t Ui + u • Vuj - .4u + V Xi p - ad Xj Eij = aE jk d Xj E ik - Y^ Au ~ K (p) d ^P> ( 3 - 2b ) 
B t + u^V£-Vu = Vu£, (3.2c) 



8 XIANPENG HU AND DEHUA WANG 

with 

K{P) '= (f+ P p)P>(l) ~ *' ^=M + (A + M)Vdiv, a = ^ ry 

We remark that for simplicity of the presentation, we will assume that a = 1 for the rest 
of this paper. 

For s £ K, we denote 

A.'f:=T-\\Z\T(f)). 

Let 

d = A _1 divu 
be the "compressible part" of the velocity, and 

u) = A _1 curlu, with curl(u))] = d X] u l — d Xi VL> 
be the "incompressible part". Setting u = A + 2/i, then system (|3.2p can be rewritten as 
pt + Ad = — pdivu — u • Vp, (3.3a) 

d t d - uAd -kp-TE = A-Miv (-u • Vu + E jk d X] E ik - j^^u - K{p)d Xi p^j , 

(3.3b) 

d t uj - pAuj -KE = A _1 curl ^-u • Vu + E jk d X] E ik - YJ^ Au ~ K (p) d ^P^j > ( 3 - 3c ) 
Ek + u • VE - Vu = VuE, (3.3d) 
u = -A -1 Vd + A^curlw, (3.3e) 

where 

1Z = A _1 curldiv, T = A~ 1 divdiv. 

The operators A, T and 1Z are differential operators of order one. 

Notice that the condition div(pF) = for all t > implies that '^j'^'J : = for all t > 
and smooth functions p, F. Hence, we have 

-1 ( d2E ij 



TE = A" 



dxidxj 
}2r 



A I . — ) - A_i divdiv(p/ + pE) (3.4) 



=o 

Ap - A^divdivQaE), 



where 



0, if i^j; 

1, if i = j. 



COMPRESSIBLE VISCOELASTIC FLOWS 9 

On the other hand, according to Lemma 16.21 (see Appendix), we have 

= A " (*> (S) - d - (§)) + ^ e -^ - 

— A 1 d Xk (Ei k ViEji — EiiViEj k ) 

= —A(Eij — Eji) + A~ 1 d Xk (EikViE i j - EijViE ik ) 

- A^d^E^ViEji - EuViEjk). 
Thus, we finally obtain 

pt + Ad = — pdivu — u • Vp, (3.5a) 

dtd - uAd - 2Ap = A- X div ^-u • Vu + E jk d Xj E ik - y^ Au ~ K (p) d ^P ~ div(p£)^ , 

(3.5b) 

d t u - pAu + A(E - E T ) = A -1 curl ^-u • Vu + E jk d Xj E ik - Y^ Au ~ K (p) d *ip) + 5 > 

(3.5c) 

(E T -E) t + u- V(E T -E)+Au = (VuE) T - VuE, (3.5d) 

St + 2Ad = -A~ x d x A- 1 ^. (u • V(Eij + E^)) + A" 1 ^ A~ l d X] ((Vu£) tJ + (VuE)^) , 

(3.5e) 

u = -A -1 Vd + A -1 curL;, (3.5f ) 

where the antisymmetric matrix S is defined as 

Sij = A 1 d Xk (Ei k ViEij — EijViE ik ) — A 1 d Xk (Ei k ViEji — EuViEj k ), 
and the scalar function S is defined as 

Sij = A l d Xi A l d X] (Eij + Eji). 
Notice that from Proposition 12. 1\ we deduce that 

||£||,b s ~ \\E + £ ,T ||ss, 

and 

\\S\\ B s + \\E - E T \\ B s m \\E\\ B s. (3.6) 
Also, according to fj3.4j) and the second equation of (|3.5j) . we have 

d t d - uAd - 2AS = A _1 div ^-u • Vu + E jk d Xj E ik - YJ^ Au ~ K (p) d ^P + div(pE)\ . 

9 (3 - 7) 
The motivation to write the second equation of (I3.5P as (I3.7P is to obtain the estimate on 

the symmetric part S of the deformation gradient, as we will see in section 4. 

The fact that this new formulation ( (|3.5h . supplemented with ()3.7p ) is equivalent to 

(|3.ip requires some explanation: it is not immediately obvious that the second and the 
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third equation in the above system are equivalent to the second equation in (|3.1|) under the 
condition that div((l + p)(I +E) T ) = 0. Notice however that the right hand side (denote it 
by O) of the third equation in (|3.5p and the term A(E T — E) are skew-symmetric matrices 
and satisfy the following Jacobi relation: 

d Xi O* + d Xj O l k + d Xk O{ = for l<i,j,k<N. 

Since ujq := A~ 1 curluo is a skew-symmetric matrix satisfying the Jacobi relation, this is 
also the case for uj. We therefore have the equivalence 

u = — A~ 1 Vd + A _1 curLj 44> divu = Ad and curlu = Au, 

which enables us to conclude that u indeed satisfies the second equation (|3.ip as soon as 
d and u satisfy the second and the third equation in (|3.5j) . 

The existence of a solution to (jl.ip is proved thanks to a classical (and tedious) iteration 
method: we define a sequence of approximate solutions of (|l.ip which solve a linear 
systems to which Proposition 14. 1 1 applies . For small enough initial data, we obtain uniform 
estimates so that we can use a compactness argument to show the convergence of such an 
approximate solution. Refer to Section 5 for more details of the complete proof. 

Let us now introduce the functional space which appears in the global existence theorem. 

Definition 3.1. For T > 0, and s £ 1, we denote 

<& 8 T = {(p, u, E) G too, T; B s+1 ' s ) n C([0, T];B s ~ l ' s ) 

x (L 1 (0,T;B 9+1 )nC([0,T\;B'- 1 )) N 

l/n t. £is+l,s\ r»/m TH. DS-l.Sl 



and 



L L (0,T;B 8+L > 8 )nC([0,T\;B a 



\\(p,U,E)\\fB«, = ||p||i«(s»-i>«) + [I u l[x,|? + \\ E \\ L ° 



NxN 



+ IIpIIl^S^+I^) + ll U llxl,(B»+l) + 11^11x1,(^+1.=-) ■ 

We use the notation QS S if T = +oo by changing the interval [0, T] into [0, oo) in the 
definition above. 

Now it is ready to state our main result: 

Theorem 3.1. There exists two positive constants 7 andT, such that, if p^ — l S B 2 '2 } 
Uo £ .B 2 1 ) f -/G-B2 !' 2 satisfy 

• ||po - l|l 5 4-i,4 + ll u oll B 4-i + ll F o - / ll i ^4-i.4 < 7; 

• div(p F([) = 0; 

• F lfc (0)V SBl F ii (0)=Fy(0)V Xj F ifc (0), 

then system (jl.ip has a solution (p,u,F) with (p — l,u, F — I) in *B~ satisfying 

\\(p- i,Q,f - 1)\\^ < r (||p - ill^^^ + l|uoll fl ^-i + I|f - 

Remark 3.1. The solution in Theorem 13.11 is also unique, but we omit the proof of the 
uniqueness. The proof of uniqueness will be same as in [3] with a slightly modification 
due to the deformation gradient. See also [20] for a proof. 
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4. Estimates of an Linear Problem 

In this section, we consider the following auxiliary linear system: 

dtp + u ■ Vp + Ad = £, 
d t d + u • Vd - uAd - 2Ap = OJl, 
< dtU + u-Vuj- /iAw + A(E - E T ) = % (4.1) 
d t {E T -E) + u- V(E T -E)+Alu = Q, 
£ t + u-V£ + 2Ad = H, 

where £, Wl, 9?, £}, 3, .ft, and u are given functions, and ()3.6p gives a relation of £7, E — E T , 
and £. We remark that, as in (13. 7p . to obtain estimates on £ , we need to rewrite the 
second equation of (I4.ip as 

d t d + u • Vd - uAd - 2A£ = 5, (4.2) 

under the constraint 

2Ap + Tl = 2A£ + 2- (4.3) 
For this system, we have the following estimate: 

Proposition 4.1. Let (p,d,to,E — E T ,£) be a solution of (|4.ip on [0,T), and 

:=jf ||u( S )|| Bf+1 d S . 
Under the condition (|4.3|) . the following estimate holds on [0, T): 

11/^)11^,4 + 11^)11^4^,4 + 11^)1^4-! + lk(t)ll B 4_ 1 

+ / (||p(s)II^4 +1 ,4 + \\d(s)\\ B N +1 + ||S(s)||^4 +li 4 + |Ks)|| B 4 +1 ) ds 



< L7e cy W{||p |||^4_ 1 ,4 + ||^o||^4_ 1 ,4 + ||do|l B 4_i + l|wo|l B £_i 
+ / e- cl " (fl) f||£||.iv 1 N+\\m\\ N t + WVtW N ,+11011-JV lN 

+ ||0:|| B 4_ 1 + p|| B 4_ 1 ,4)d s }, 

where C depends only on N . 

Remark 4.1. Notice that the constraint (|4.3p is always satisfied by our system (|3.5p in 
view of the divergence property of the deformation gradient E, and 5 will be given as the 
right-hand side of ()3.7p . This implies that the estimates in Proposition 14.11 also hold for 
solutions to ()3.5p . 

Proof of Proposition \4-l\ To prove this proposition, we first localize (|4.ip in low and high 
frequencies according to the Littlewood-Paley decomposition. We then use an energy 
method to estimate each dyadic block. To this end, we will divide our proof into four 
steps. 

Let (p,d,u,E) be a solution of (I4.ip and K > 0. Define 

p = e- KV ®p, d = e- KV ^d, E T -E = e- Kv ^{E r - E), 
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and 

Applying the operator A q to (]4.ip and (14. 2j) . we deduce that (A g p, A 9 d, A g w, A q E) satisfies 

' 8 t A q p + A 9 (u • Vp) + AA 3 J = A 9 £ - ^'(t) A,p, 
d t A q d + A g (u • Vd) - vAA q d- 2AA q p = A q Wl - KV'(t)A q d, 
d t A q u + A q (u ■ Vu) - jjAA q u + AA q (E - E T ) = A q m - KV'(t)A q u}, 
d t A q (E T -E) + A g (u • V(E T - E)) + Aoj = A q Q - KV'(t)A q (E T - E), 
d t A q d + A q (u ■ Vd) - uAA q d- 2AA q £ = A q Z, 
k d t A q £ + A 3 (u ■ V£) + 2AA q d = A q k. 

Denote 

9q : = 2«(t -i) (2||A ?/ 5||£ 2 + 2||A 9 d~||£ 2 + ||A 9 (£ T - E)f L2 + ||A g £||| 2 + \\A q oj\\ 2 L2 



- —(AA q (E - E)\A q u>) - ~(AA q p\A q d) - -(AA q £\A q d)) 

for q < q with rj = max{^, u, ^} + 1; 

g q := 2«(t -i) (||AA 9 p||| 2 + ||AA g (E T - £)||| 2 + ||AA 9 £||| 2 + ||A g J||| 2 + ||A g cD||| 2 
- (AA ?y 5|A ? d) - (AA 9 (£ T - S)|A,w) - (AA,£|A,d)^ 

for q > go, where /3i = /?2 = -jr, 7 = max ^ -7 > + 1, and go is chosen to satisfy 

||AA,/|| L 2 > 2 7 ||A f/ /|| L2 for all q > q . (4.5) 
Due to the fact suppJ r (A g p) C 2 q C and suppJ r (A q E) C 2 9 C, one deduce that 



< C (4.6) 

V2^- 1 ^W(||A g p|| i2 + \\A q E\\ L .) + 2^T-D(||A,d|| ia + ||A^|| i2 )7 

for a universal constant C. 

The first two steps of the proof are devoted to getting the following inequality: 

\±- t g\ + k2^ N ' N ~ 2 ^ (\\A q p\\\ 2 + ||AA 9 (£ T - E)\\\ 2 + \\AA q £f L2 + ||A^||| 2 + || A g d~||| 2 ) 

< c aq9q (\\i\\^ _ lif + H0ii 5f _ ljf + \\k\\ § n + \m B N^ + m\ B ^ + \m B %-i 

+ V'(t)(\\p\\ B *_ l4 + \\E T - £|| -^,4 + 11^11^^,4 + IMI^ + IHI B £-i)) 
-KV'gl 

(4.7) 

where k is a universal constant. 
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First Step: Low Frequencies. Suppose q < q$ and define 

f 2 q = 2\\A q p\\ 2 L2 + 2||A 9 d||| 2 + \\A q (E T - E)f L2 + ||A,£ 1|| 2 + ||A ? u,||| 2 

- — (AA q (E T - E)\A q Cj) - -(AA q p\A q d) - -(AA q £\A q d). 
rj 77 77 

Taking the L 2 -scalar product of the first equation of (|4.4p with A q p, of the second equation 
with A q d, of the third equation with A q cj, and of the fourth equation with A q (E T — E), 
we obtain the following four identities: 

— IIA^Hia + (A 9 (u • Vp)\A q p) + (AA q d\A q p) = (A q £\A q p) - KV'\\A q p\\ 2 L2 ; (4.8) 

— llMlla + ^l|AA ? d||| 2 + (A q (u • Vd)\A q d) - 2(AA q p\A q d) 
= (A q m\A q d)-KV'\\A q d\\ 2 L2 ; 

^l|A^|| L2 + (A 9 (u • VQ)\AgU>) + p\\AA q u\\ 2 L2 + (AA q (E - E T )\A q ti) 
= (A q m\A q Co) - KV'\\A q u\\ 2 L2 ; 

lj t \W(E T - E)\\ 2 L2 + (A 3 ((u • V(E T - E)))\A q (E T - E)) + (AA q ub\A q (E T - E)) 
= (A q (Q T - Q)\A q (E T - E)) - KV'\\(E T - E)\\ 2 L2 . 

(4.11) 

And, we also have, from fifth and sixth equation in (|4.4p 

lj t \\A q df L2 + i/||AA 9 d||| a + (A,(u • Vd)\A q d) - 2(AA q £\A q d) 
= (A q 3\A q d) - KV'\\A q d\\ 2 L2 ; 

\j t \\\£\\h + ( A «( u • V£ )|A/) + 2(AA q d\A q £) = (A q k\A q £) - KV'\\£\\\ 2 . (4.13) 

For estimates of the term (AA 9 p| A q d), we apply A to the first equation in (|4.4p and 
take the L 2 -scalar product with A q d, then take the scalar product of the second equation 
with AA q p and sum both equalities to yield 

j t {AA q p\A q d) + HAAgdll^, - 2||AA f?/ 5|| 2 2 + (AA q p\A q (w ■ Vd)) 

+ (AA,(u • Vp)\A q d) + u{A 2 A q d\AA q p) ( 4 - 14 ) 
= (AA q £\A q d) + (AA q p\A q m) - 2KV'(AA q p\A q d). 

For estimates of the term (AA q (E T — E)\A q u), we apply A to the fourth equation in 
(|4,4p and take the L 2 -scalar product with A q u), then take the scalar product of the third 
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equation with AA C/ (£' T — E) and sum both equalities to yield 

^(AA 9 (£ T - E)\A q u) - \\AA q {E T - E)\\\ 2 + ||AA g £||| 2 + (AA q (E T - E)\A q (u • Vw)) 

+ (AA,(u • V(£ T - S))|A,w) + /z(A 2 A g LD|AA g (£ T - E)) 
= (AA q O.\A q Co) + (AA q (E T - E)\A q W) - 2KV'{AA q {E T - E)\A q u>). 

(4.15) 

For estimates of the term (AA g £| A q d), we apply A to the last equation in (|4.4|) and 
take the L 2 -scalar product with A q d, then take the scalar product of the fifth equation 
with AA q £ and sum both equalities to yield 

j t (AA q £\A q d) + ||AA 3 d|| 2 2 - 2||AA f? £|| 2 2 + (AA 9 £|A g (u • Vd)) 

+ (AA g (u • V£)\A q d) + u(A 2 A q d\AA q £) ( 4 - 16 ) 
= (AA q k\A q d) + (AA q £\A q 3) - 2KV'(AA q £\A q d). 

Taking linear combination of (|4.8p - (|4.16p . we obtain, 

Ydft + ^( 2 II AA ^I" 2 + H AA «^ T " V)t L , + 2||AA/||| 2 + (2 V - 3)||AA 9 d|| i2 
- (i>A 2 A q d\AA q p) - (pA 2 A q u\AA q (E T — E) — i>(A 2 A q d\AA q £ ))) + KV 'f q 



= x, 

(4.17) 

where 

X := 2(A q £\A q p) + (A q M\A q d) + (A q m\A q u) + (A q Q\A q (E T - E)) + {A q Z\A q d) 
+ (A g £|A 9 £) - 2(A g (u • Vp)\A q p) - 2(A g (u • Vd)\A q d) 

- (A 9 (u • V(£ T - -E))|A g (£ ,T - E)) - (A g (u • V£)|A g a>) - (A g (u • V£ )|A 3 £) 
+ ^[{AA q p\A q (u ■ Vd)) + (AA,(u • Vp)\A q d) + (AA 9 (£ T - E)\A q (u ■ V£)) 

+ (AA g (u • V(£ T - £))|A,£>) - (AA g £|A g d) - (AA 9 (£ T - E)\A q % 

- {AA q p\A q m) - (AA q Q\A q Cj) + (AA g £|A 9 (u • Vd)) + (AA g (u • V£)\A q d) 

- (AA q K\A q d) - (AA q £\A q J)Y 

As q < qo, there exists a constant Co > 1, which depends on X, p, such that 

I/ 2 < \\A q p\\h + \\A q d\\l 2 + \\A q (E T - E)\\ 2 L 2 + ||A^||| 2 + ||A^||| 2 < c / 2 , (4.18) 
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and this equivalence implies that there exists a universal positive constant k depending 
on A and fx, such that 

^{2\\AA q p\\ 2 L2 + \\AA q (E T - E)\\ 2 L2 + 2||AA (? £|| 2 2 + (2r? - 3)||AA g <2|| L 2 



2'/ 



{uA 2 A q d\AA q p) - (pA 2 A q Cj\AA q (E T - E) — u(A 2 A q d\AA q £)) 



(4.19) 



> K2 2g (j|A g/ 5|| 2 2 + ||A ? (£ T - E)f L2 + \\A q £\\ 2 L2 + ||A,w|||a + ||A,d||] 

For terms on the right-hand side of (|4.17p . we use Lemma 12.21 (|4.18p . and the Cauchy- 
Schwarz inequality to obtain 

\X\ < (7/ ? (||A ? £|| L 2 + ||A,07t|| L 2 + ||A g «tt|| L 2 + ||A ? 0|| i2 + \\A q k\\& + \\A q 3\\ L 2 

+ 2-^-Va q V'(\\p\\ § N_ li N + \\E T - E\\ 6 »-i,» + \\£\\ s %-x,% ( 4 - 2 °) 

+ NlU_ 1 + Nl fl4F _ 1 )). 

Hence, combining (liTfl) . KIM . KW\ . (1420)) together, we obtain 

+ Kl 2 ^\A q ~pf L2 + ||AA 9 (i T - E)\\\ 2 + \\AA q £\\\ 2 + \\A q u\\ 2 L2 + ||A ? d||| 2 
< Cf q (\\A q Z\\ L2 + ||A 9 Stt|| L 2 + IIA^H^ + ||A,0|| L 2 + ||A^|| L 2 + \\A q 3\\ L 2 



(4.21) 



+ 2-^-^a q V'(\\p\\ § M^N + \\E T - E\\ § *_ h N + \\£\\ § N_ h N + \\d\\ B N_ 
+ \\u\\ B ,_ 1 ))-KV'f 2 . 



Second Step: High Frequencies. In this step, we assume q > qo. We apply the 
operator A to the first equation of (|4.4|) . multiply by AA q p and integrate over R N to yield 

lj t W AA iP\\ 2 L* + ( AA <?( U • v P)\ AA iP) + (A 2 A,<l|AA 9 p) ^ 
= (AA q £\AA q p) - KV'\\AA q p\\ 2 L2 . 

Applying the operator A to the fourth equation of (|4.4|) . multiplying by AA g (.E T — E) 
and integrating over R , we get 

^||AA g (£ T - E)f L2 + (AA g (u • V(E T - E))\AA q (E T - E)) + (A 2 A,u;|AA g (£ T - E)) 

= {AA q a\AA q {E T - E)) - KV'WAA^E 1 - E)\\ 2 L2 . 

(4.23) 
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Applying the operator A to the sixth equation of (j4.4|) . multiplying by AA q £ and inte- 
grating over WL N , we obtain 

lj t \\AA q Sf L2 + (AA q (u ■ V£)\AA q £) + 2(A 2 A q d\AA q £) 
= (AA q k\AA q £) - KV'\\AA q £\\ 2 L2 . 

Denoting 

f q = \\AA q ~pf L2 + \\AA q (E T - E)\\ 2 L2 + ||AA g 5||| 2 + 2 7 ||A ? d||i 2 + 7 ||A,a;||| 2 
- Pi(AA q p\A q d) - j3 2 (AA q {E T - E)\A q £j) - 2^{AA q £\A q d). 

Combining (gSJ, (ITOD . KT2\\ . (gUP , (fl~T5l) . g22|), and (IMjl . we obtain 

^/ 9 2 + /3 1 ||AA 9 p||| 2 +2/3 1 ||AA^||| 2 + (2^-^) l|AA ? d||| 2 + (j*7 - y) l|AA g u)||| 2 

+ Y'||AA g (E ,T - £)||| 2 - 2 7 (AA 9 p| A 9 d) - 2 7 (AA^|A g d) - 7 (AA f/ (£ T - E)\A q co) 
+ KV'f 2 q 

= y, 

(4.25) 

where 

y = i(A q m\A q d) + 7 (A q 3\A q d) + 7 (A g 0T| A 9 cD) - 2 7 (A,(u • Vd)\A q d) 

- 7 (A,(u • Vu)\A q uj) + ^ ((AA,p| A,(u • Vd)) + (AA 9 (u • Vp)\A q d) - (AA q £\A q d) 

- {AA q p\A q Wl) + 2(AAg£|Aq(u • Vd)) + 2(AA,(u • V£)\A q d) - 2(AA q k\A q d) 

- 2(AA q £\ A 5 5))) + y ((AA 3 (£ T - £)|A g (u • Vw)) + (AA g (u • V(E T - E))\A q u) 

- (AA q Q\A q Cj) - (AA q (E T - E)\A q Vlj) - (AA q (u ■ Vp)\AA q p) 

+ (AA q £\AA q p) - (AA q (u ■ V(E T - E))\AA q (E T - E)) + {AA q Q\AA q (E T - E)) 

- (AA 9 (u • V£ )\AA q S) + (AA q k\AA q £). 

Notice that, for q > qo, we have 

f q « ||AA 9 p||| 2 + ||AA g £||| 2 + ||A 9 d|| 2 2 + ||A^|| 2 2 ; (4.26) 

and 

Al|AA 9 p||| 2 + 2/3 1 ||AA/||| 2 + (Try,- ^ ||AA ? d||| 2 + (m7~ y) l|AA 9 o>||£ 2 

+ ^|||AA ? (£ T - £)||| 2 - 2 7 (AA 9 p|A,d) - 2 7 (AA ? £|A ? d) - r y(AA q (E T - E)\A q O) 
a ||AA^||| 2 + ||A g d||| 2 + ||A^||| 2 + ||AA 9 (£ T - E)\\l + ||AA/||| 2 . 

(4.27) 
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Next, we apply Lemma 12.21 to obtain, using (|4.26p 

\y\ < Cf q (\\AA q 2\\ L * + ||AA 3 0|| i2 + ||AA^|| La + \\A q 3\\v + \\A q M\\ L2 + ||A 9 <tt|| L2 
+ a q 2-<*f -QVQIpW -^.^ + lldl^., + \\E T - EW^n 

+ \\i\\ § N_ 1 N +|H| B 3_ 1 )) . 

(4.28) 

Therefore, there exists a universal positive constant k, which depends on /i and u, such 
that, 

i|/ 9 2 + K (||AA 9/ 5||| 2 + ||AA^||i 2 + ||A 9 d||| 2 + ||A ? a}|| 2 2 ) 

< Cf q (\\AA q £\\ L2 + ||AA g £)|| L2 + \\AA q k\\ L 2 + \\A q 3\\ L 2 + ||A,SOT|| L 2 + \\A q m\\ L 2 
+ a^^-Vv'm^n + ||J|| b4 _ 1 + \\E T - E\\^_ h N 
+ \\S\\^ 1 ,+\M B ,_ 1 ))-KV'fl 

(4.29) 

Third Step: Damping Effect. We are now going to show that inequality (|4.7p 
entails a decay for p, E, d and uj. We postpone the proof of smoothing properties for d 
and uj to the next step. Let 5 > be a small parameter (which will tend to 0) and denote 
h q = 9q + <5 2 - From (|4.7p . and dividing by h q , we obtain 

d u 

-h q + Kh q 

< Ca q (\\£\\ § N_^n + .^^^ + u + WZW^ + nscrtii^ + \\%\ B %-i) 

+ Ca q V'(\\p\\^ + || d|| B 4 + ||£ T - ^l^* + ||£|| B g: + ||^H Bf _ x ) 

- Ary'/ig + 5Ary' + <fc. 

Integrating over [0, t] and having 5 tend to 0, we obtain 

9q(t) + « / 9q(r)dT 
Jo 

< g q (0) +C [ a q ( S )(\\Q\\-j* + P||^_ lif + 1151^^ + ll^ll B f -x + IWI B f -x 

+ ||0|| Bf + jf V'( S ){Ca q (s)(\\p\\^ _ ljf + Hdll^., + ||£ T - ^H Bf _ lif 

+ Hfll^.!,^ + \\ti\\ B N_^ - Kg q (s)}ds. 

(4.30) 
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Thanks to (|4.6p . we have 



Ca q (s)[\\p\\^N_^N + ||d|| B £-i + \\E - E\\.N_ h N + 11^11^4-1,^ + 11^11^4-1 J - K 9g{s) 



< q (s)\\E T (s)-E(s)\\„»_ lS -i 

N -i N 



«_, « 

+ Ca q ( S )\\E T (s)-E(s)\\ i .,_ 1 ,-KC- 1 2^- ,T ^\\A q (E T - E)\\ L * 



+ Ca q (s)\\d(s)\\ D ,_ 1 -KC- 1 2^- 1 ^A q d\\ L 



B"2~ 



2 



+ Ca q {s)\\u:{ S )\\ B N_ 1 -KC-^- 1 ^\A q u\\ L2 . 



If we choose K > C , we have 



J2 {Ca q (s)(\\p\\^N_ lt N + \\d\\ B %-i + ||£ T - E\\^_ h N + ||£|U 



~ iV 1 iV 



11 "BT 



-^<7 9 ( S )} <0. 



According to the last inequality, and thanks to (|4.6p and (|4.30|) . we conclude, after 
summation on Z, that 

||/3(i)|| sf + + ||^(i)ll B f -i + \\d(t)\\ B f-i 

< c{\\ P0 \\-N_ h N + \\E \\-N_ lt N + ||wo|| B «: _! + \\M B %-i + / (l|£(s)ll B £_i,f 

+ \\m(s)\\ B N , + ||5(s)|| Bf _! + PII B f - 1: f + l|9fa)ll fl S-i + ll^ll B f )<**}■ 

(4.31) 

Fourth Step: Smoothing Effect. Once showing the damping effect for p and E, 
we can further get the smoothing effect on d and oj for the system (|4. 1[) by considering 
the second and third equations with terms Ap and AE being seen as given source terms. 
Indeed, thanks to (I4.3ip . it suffices to state the proof for high frequencies only. We 
therefore assume in this step q > qo- Define 

I q = 2^f- 1 )||A,(i|| L2 +2^ 1 )||A g cD|| L2 . 
Then, from the energy estimates for the system 

jd t A q d + A q (u ■ Vd) - vAA q d = 2AA q p + A q Wl - KV'(t)A q d, 
\ d t A q uj + A 9 (u • VlD) - pAA q u = TZA q E + A q m - KV'(t)A q u, 
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we have 

- — T 2 4- K2 2q T 2 

< I q M\\A q p\\v +2«(f- 1 )||A^|| L2 +2^-V\\A q m\\ L 2 + 2"f \\A q E\\ L2 



+ I q V\t)(Ca q {\\d\\ B N_ 1 + \\u\\ b n_ x ) - KI q ), 



for a universal positive constant k. Using J 2 = I 2 + 5 2 , integrating over [0, t] and then 
taking the limit as 5 — > 0, we deduce 



I q (t) + K2 2q [ I q (s)ds < 7,(0) + f 2 <1 ^- 1 \\\A q Wl(s)\\ L 2 + \\A q m(s)\\ L 2)ds 
Jo Jo 

+ f 2^{\\A q p{s)\\ L 2 + \\A q E{s)\\ L 2)ds (4.32) 
Jo 

J 

We therefore get 

Y, 2^- l \\\A q d{t)\\ V2 + \\A q uj{t)\\ L .) + « f E 2 9( ^ +1) (l|A g (i( S )|| L2 + ||A^(t)|| L2 )d S 
g>go g>go 

<\\do\\ B N_ 1 + \\u \\ B N_ i + f (Wmis^^^ + lMs^l^^ds 

J 

+ / E2^(||A ? p( S )|| i2 + ||A g 75( S )|| I , 2 )^ + Cn*) sup(||d|| ^_ 1 + |H| * J. 
Using (|4.3ip . we eventually conclude that 

K f Y 2^ +1 \\\A q d(s)\\ L 2 + \\A q u(s)\\ L2 )ds 

< C(l + U(t)){||p ||^4- 1 ,4 + ll^oll^-^ + llwoll^.! + ||do|| fl 4_i 

+ f (\\Z(s)\\-N_ h N + \\M(S)\\ B %-X + \M^)\\ B ^-1 + I^WII^-l^ 
J 

+ \\&(s)\\~N_ h N + \\3( S )\\ B N_ 1 y s y 

Combining the last inequality with (|4.3ip . we finish the proof of Proposition 14.11 □ 

In the rest of this section, we will sketch the proof of the existence of a unique global 
solution to (|4.ip . To this end, we only need to use some properties of transport and heat 
equations in nonhomogeneous Sobolev spaces TP with high regularity order, see [TJ [2]. 
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Indeed, we set (p° ,d° ,uj° , E°) = (po, do, ujo, Eq) and 

' d t p n+l + u • Vp n+1 = -Ad n + £, 
d t d n+l - isAd n+1 = Wl - u • Vd n + 2Ap™, 
< d t 0J n+l - pAu n+1 = m + u • Vw" + A((£ n ) T - £ n ), (4.33) 
d t ((E n+1 ) T - E n+1 ) + u • V((E n+1 ) T - E n+l ) = -Auj n + 0, 
d t £ n+1 + u • V£ n+1 = -2AcT + £, 

where 

{p n+ \d n+ \u n+ \E n+l )\ t=G = (p ,d ,u ,E ) 

Let T > and s (large enough) be fixed, and let K be a suitably large positive constant 
(depending on s, T and u). We set 

p n = e - Kt p n , d n = e~ Kt d n , co n = e- Kt u n , E n = e~ Kt E n . 

The straightforward computations show that {(p n , d n , to 11 , E n )} n£ ^ is a Cauchy sequence 
in 

C([0,T];H s ) x (CQO.T];/^ 1 x H 3 " 1 ) n L\[0,T]; H s+1 x H S+1 )) 2N x C([0,T]; H s ) NxN . 

Denoting by (p,d,Cb,E) the limit, then it is easy to show that (e Kt p,e Kt d,e Kt Cu,e Kt E) 
solves (|4"TTj) . 

5. Global Existence 

The goal of this section is to prove the global existence of solutions to fjl .3j) by build- 
ing approximating solutions (p n ,u n ,E n ) using an iterative method. Those approximate 
solutions are solutions of auxiliary systems of (|4.ip to which Proposition 14.11 applies. 

We set the first term (p°,u°,E°) to (0,0,0). We then define {{p n , u n , E n )} nGM by 
induction. In fact, we choose (p n+1 , u n+1 , E n+1 ) as the solution of the following system: 

'd t p n+l + u n • Vp n+1 + Ad n+l = £ n , 
d t d n+l + u n ■ Vd n+1 - vAd n+1 - 2Ap n+l = Wl n , 
d t oj n+1 + u" • Vuj n+l - pAuj n+1 - A(E n+1 - (E n+1 ) T ) = m n , 
< d t ((E n+1 ) T - E n+l ) + u • V({E n+1 ) T - E n+l ) + Auj n+1 = Q n , (5.1) 
d t £ n+1 + u • V£ n+1 + 2Ad n+1 = & n , 

u n+l = _ A -l vd n+l _ A- 1 curlw n+1 , 

Xp n+ \d n+ \uj n+ \E n+1 )\ t=0 = (^,A- 1 divu„,A-Wu B ,^). 

where 

Pn = ^ ^<lP0i u n = ^2 A 9 U 0' E n = \ E 0, 
\q\<n \q\<n |<?|<n 

= -p n divu n , 

m n = u n • V<F - A-Mivf u n • Vu n + K(p n )Vp n + — " „4u" 

V 1 + p n 

-E? k d Xj E? k + div(p n E n )), 
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= u n • Vw" - A-Wl ( u 11 ■ Vu n + K(p n )Vp n + — ^Au n - E%d x El) , 

\ 1 + p U J 3 J 

Q n = (Vu n E n ) T - Vu n E n , 

and 

= u" • V£ n - A^A" 1 ^. (u • V(El) + 

+ A~ 1 Sl Bj A- 1 flL j ((Vu B ^ + (Vu'T)^). 

As in (|3.7p . due to the divergence property of the deformation gradient, we can rewritten 
the second equation in f)5. 1 1) as 

5 t d n+1 + u • VcT +1 - uAd n+1 - 2A£ n+1 = T, (5.2) 

where Z n is given by 

T = u n - Vd n - A-Mivfu" • Vu n + K(p n )Vp n + -Ji—Au n - Ej k d Xj E? k - div(p n E n 

The argument in the previous section guarantees that the system (|5.ip is solvable and 
Remark 14.11 tells us that, in view of (|5.2p . the solution to (|5.ip satisfies the estimates in 
Proposition 14.11 

5.1. Uniform Estimates in the Critical Regularity Case. In this subsection, we 
establish uniform estimates in 53 ~ for (p n ,u n ,E n ). Denote 

7 = Po ~ i jv + Un jv 1 + En - jv i jv . 

We are going to show that there exists a positive constant T such that, if 7 is small enough, 
the following estimate holds for all n G N: 

11(^,^,^)11^, <r 7 . (*p„) 

We will prove (^3 n ) by the mathematical induction. Suppose that (^3 n ) is satisfied and let 
us prove that also holds. 

According to Proposition 14. II and the definition of (p n , u n , E n ), the following inequality 
holds 

)|( re+1 n+1 n+1) H N <Ce cvn (\\p \\-N ^H-Huoll jv , + IISqII-jv X iv 

+ ||£™|| , -JV 1 N+\\Wl n \\ JV lN + ||ari| , JV u 

(5.3) 

where 

poo 

v n = y o iiu"( S )n sf+1 d S . 

Hence, it remains to obtain the estimates for £ n , 9Jt n , 9T n , & n , Z n and n by using (*$«). 

Estimate of £ n : The estimate of £ n is straightforward; thanks to Proposition 12.21 we 
have 

||£ n || ~ jv lN <C||p n || -jv lN J|divu n || t/ « 

< cry. 
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Estimates of 5D? n , 9T n , and 2 n : To estimate 9JT n , OT" - , and Z n , we assume that 7 
satisfies 

1 

7< 



2re 



where £ is the continuity modulus of B 



iv ^ JV 



and 



11 ^11 ^ (f\\ „ n \ 

IIP Hl^r+xr^) S <-\\p I 

l-E" 1 !! T°o(W>+ vTCPAM < C||-E n 



L°°. If (<p n ) holds, then 
1 

JV 1 Jv < — , 



1 

< -. 



TR+xR-N") - -11- n B ^—i,^ — 2 

We will estimate 9K n term by term. First we have, according to Proposition 
Lemma 12, 11 



1 + p" 



-V 2 u n 



JV , 



< C||V 2 u n | 



P" 



l + p n 



L°°(BT) 



and 



(5.5) 



2„,2 



< C||u n || , N+1 \\p n \\ n < CT z j 



Since K(0) = and Lemma 12.1^ one has 



„n\\2 



\\K(p n )Vp n \\ , N ! < C||ir(p n )|| „ ^J|V/) n || „ JV ! <C||/|| 

II L'fBT" ) - 11 VF ;ll L 2 (_B^-) n r M L 2 (BT _1 ) _ 111 ll L 2 (B%) 

On the other hand, we have, by Holder's inequality, 

. _ / JV, JV , , \ 2 / JV,, JV , , 

£ 2^" ,T W||A s p»(*)|| £a ) 2*^ + "<*>||A,p» 



P JV 

IK V(B^) 



dt 



< 



~ JV i JV P 



B^" + i '"Z" 



I 1 1 ll„?i|| ^ /^"p2.,2 

\P \\ - Jv , jv p , - jv ,, jv <01 7 . 
lr L°°(B"2" _ '"2") "Li(B"2" +1 '"2") ~~ 



I jv , < cr 2 7 2 . 



Thus, the above two inequalities imply that 

\\K(p n )V P n 

From Proposition 12.21 we obtain the following estimates: 

\\u n -Vd n \\ , jv , +||u n -Vu"|| jv , +||u n -Vw n || n , 

Icf 



(5.6) 



< C 

< C||u n | 



u 



B"2"~ 



JV + \\L0 



jv , , dt 



L°°(B~2~ ) 



i \\d n \\ , jv +1 + ||w n || , 



< CT 2 7 2 ; 



\E 1 Aud x .E 1 ! i h \\ jv , < C 

i j/e J WCM Z, 1 (B"2" ) ~~ 



IE 1 '™ 



L 1 (B^" +i ) 



(5.7) 



B~2~ '"2~ 



- JV,, JV (it 
BT +1 T 



^ /"Ml IT" 1 II II T? n II 

^ O Lfi - jv , jv, jo -,,-JV,iJv, 

— 11 "loo(BT _1 'T) ll L 1 (B^" + 1 '^") 

< CT 2 7 2 ; 



(5.8 
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and 



oo 







c(llp n ll- 



- N , N \\£L , ~ N ,1 N 
Loo(S"2" _1 '^")" N L1(BT +1 'T) 



_L 1 1 IT" 1 II 1 1 rP' 1 1 

+ Lc/ -AT . AT lO , - N N 



Hence 



\Q n \\.N ! iV + -iV , iV < CT 7 . 



From ([O]) , (pU0|) and ((IDT]) , we finally have 



Now we choose T = 4C and choose 7 such that 



(5.9) 



< ctV- 

Summarizing (|5.5p - (|5.9|) . we finally get 

\\w n \\ * u + ||w n || ,, iv ls + 1131 n u <cr 2 7 2 . (510) 

Estimates of n and ^ n : It is easy to obtain, using Proposition 12.21 

/•oo 

||Vu n £ n || -n , N <C \\u n (t)\\ N +1 \\E n (t)\\-N lN dt 

<C||u n || „ N+1 \\E n \\ -n (5-H) 

— 11 ll L 1 (_B^"+ 1 )" L°°(.B"2" '"2") 

< cr 2 7 2 . 



IKp^u^ 1 ,^ 1 )!!^ < Ce cr ^(rV + 7)- (5.12) 



r 2 7 <l, e cr ^<2 and r 7 < -L (£) 



then (<p n ) holds for all n G N. 



5.2. Existence of a solution. In this subsection, we show that, up to a subsequence, 
the sequence (p n ,u n ,E n ) n& ^ converges in V'(R + x R N ) to a solution (p,u,E) of (TQD . 
We will use some compactness arguments. The starting point is to show that the first- 
order time derivative of (p n ,u n ,E n ) is uniformly bounded in appropriate spaces. This 
enables us to apply the Ascoli-Arzela theorem and get the existence of a limit (p, u, E) 
for a subsequence. Then, the uniform bounds of the previous subsection provide us with 
some additional regularity and convergence properties so that we may pass to the limit in 
the system (|5.1|) . 

To begin with, we have to prove uniform bounds in the suitable functional spaces and 
the convergence of {(p n , u n , E n )} ne fq. The uniform bounds is summarized in the following 
lemma. 
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Lemma 5.1. {(/?", u", E n )} n ^ is uniformly bounded in 

N / i , . \ NxN 



1 / v .\ / I 



JV 3 



CL i^J x B~- ) j x (Cf oc 

i / AT 1 f N o\ N / N ,\NXN\ 

(and also in C t 2 oc I B^" 1 x f^ T J x (Ba" M J ifN>3). 

Proof. We will finish the proof via five steps. 

Step 1: Uniform bound of dtp n in L 2 (S~ _1 ). In fact, notice that 

d tP n+1 = -p"divu n - u n • Vp n+1 - Ad". 

According to the estimates in the previous subsection and the interpolation result in 
Proposition 12. 1\ we see that {u"} ne N is uniformly bounded in L 2 (B~), and {(p n , E n )} ne ^ 
is uniformly bounded in (L°°(B~)) n2+1 . Thus, — p"divu n — u n • V p n+1 — Ad" is uniformly 
bounded in L 2 (B~~ l ), which implies that dtp n is uniformly bounded in L 2 (B~~ l ), and 

furthermore {/9"} n , S N is uniformly bounded in Cf oc B~~^j. 

Step 2: Uniform bound of dtE n in L 2 {B~~ l ). In fact, notice that 

d t {(E n+1 ) T - E n+1 ) = -Ku n + Vu n E n - u" • V((£" +1 ) T - E n+l ). 
According to the estimates in the previous subsection and the interpolation result in Propo- 
sition I2TT1 we conclude that {u n } n6 N is uniformly bounded in L 2 (B^), and {(p n , E n )} n£ fq 
is uniformly bounded in (L 00 (Bt)) n2 + 1 . Thus, -Aw" + Vu"£"-u"-V((£" +1 ) T -£" +1 ) 
is uniformly bounded in L 2 (i?"2" _1 ). Therefore, dtE n is uniformly bounded in L 2 (B~~ l ), 

NxN 



and furthermore, {(E n ) T — E n } n ^ is uniformly bounded in ( Cf oc (lR + ;.B2 



1 



Similarly, we can show that {£ n } n ^ is uniformly bounded in Cf (R + ;i? 2 1 j. Hence, 



JV 

loc 

N 



{E n }nem is uniformly bounded in C^ c B 2 1 

Step 3: Uniform bound of dtd n in L5 ) + L 4 (f?~2~~i). To this end, we recall 
that 

8 t d n+1 = -u n ■ Vd n+1 + vAd n+1 + 2Ap n+1 + u" • Vtf" 

- A^div^u" • Vu" + K{p n )Vp n + Y^Avr - EJ k 8 Xj E? k + div(p n E n f) . 
The estimates in the previous subsection and the interpolation result in Proposition 12.11 

JV 1 4 JV+1 

yield that {u"} ne ^ is uniformly bounded in L°°{B^~ )nZ>3 (B^~ ). This uniform bound, 

JV 

combining together with the uniform bound of {p n } n &N in L°°(B~), gives an uniform 
bound of 

-u" • Vd n+1 + uAd n+1 + u" • VdP - A _1 div ( u" • Vu" + — ^Au n 



l + p n 

in (B^2~ ). Using the uniform bound of {p n } n eN in L°°(B~) n L 2 (B~) obtained from 

~ JV . JV -1 ~ JV 11 JV 

the uniform bounds of {/?"} n eN hi L°°{B~~ '~)C\L (B~ + '~) and the interpolation result 
in Proposition 12. 1\ we deduce that {p"}neN is uniformly bounded in L 4 (B 2 ) 5 and hence 
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JV-3 



{Ap n } n& ^ is uniformly bounded in L (B 2 ). Finally, both E^ k d Xj E™ k and div (p n E n ) are 



JV-3 



also uniformly bounded in L (B 2 ), since {p n } n ^ and {E n } n ^ are uniformly bounded 

AT . JV — 1 

in L°°(B~) n L (B~2~) and div(p.E) can be rewritten as the sum of pdwE and X7pE. 



Therefore, {SfCpj-ngN is uniformly bounded in L3 (£? 2 2) + L (B 2 



2 



Step 4: Uniform bound of c^w" in L 2 ) + L 4 (i? ^2 3 ). We recall that 
d t u n+1 = -u n ■ Vuj n+1 + JlAuj n+1 + A({E n+1 ) T - E n+l ) + u n • Vw" 



A-Wl (u n • Vu" + K(p")Vp" + j^Au n - E] k d x .E? k 



Similar to the argument in Step 3, we conclude that 



-u n • Vd n+1 + vAd n+1 + u n • Vd n - A' 1 carl ( u n • Vu n + — P - — ^u r 

V i + p n 

4 JV— 3 JV 

is uniformly bounded in Ls{B^~), using the uniform bound of {p n } n eN i n L°°(B^) n 
L 2 (B^~). Also, similarly to Step 3, {Ai? n+1 } neN and {£ , ™ fc 5 :rj -E™ fc } neN are uniformly 

bounded in L A (B^^). Hence, {dtUJ n } n eN is uniformly bounded in L% (B^~i) +L A {B^~i). 
Step 5: Uniform bound as N > 3. Indeed, in this case, the only difference is the 

JV -1 

uniform bound on u n . Actually, from the uniform bounds on {u n } ng pj in L°°(B^~ ) D 
L 2 (£>~2~), we deduce that {u n - Vu n } ng N is uniformly bounded in L 2 (B^~~ 2 ). Then, follow- 
ing the same argument in Step 3 and Step 4, we can deduce that {dtd n } n ^ and {dtco n } n ^ 
are uniformly bounded in L 2 (£>" 2 ) + L°°(i?" 2 ) ) because {Ap n } ne ^, {AE n+1 } ne ^ and 

JV o 

{Ej k d Xj E™ k } ne f$ are uniformly bounded in L°°(B 2 z ). This further implies that {5 t u n } nG N 
is uniformly bounded in L 2 (B~2~ 2 ) + L°°{B^~ 2 ), and hence is uniformly bounded in 

J iocy 



We can now turn to the proof of the existence in Theorem I3.1[ using a compactness 
argument. To this end, we denote {XpjpeN be a sequence of Cq'(U. n ) cut-off functions sup- 
ported in the ball B(0,p+1) in M. N and equal to 1 in a neighborhood of the ball B(0,p). For 
any p G N, Lemma [541 tells us that {{XpP n -> Xp u? \ XpE n )} n eN is uniformly equicontinuous 
in CXR+jfif- 1 x (B^) 1 * x (B^~ 1 ) NxN ). Notice that the operator / h-> x P f is compact 

JV_i JV ■ JV_i iV 1 JV-3 ■ JV-3 

from B 2 1 ni?2 into H 2 , and from i? 2 L nB 2 into H 2 . This can be proved eas- 
ily by noticing that / 1— > Xpf is compact from H s n H s into .fP for s < s' and -B s =— )■ 
We now apply the Ascoli-Arzela theorem to the sequence {{x p P n -, Xp u "> XpE n )} n eN on 
the time interval [0,p\. We then use Cantor's diagonal process. This finally provides us 
with a distribution (p,u,E) belonging to C(R + ;H^~ 1 x (H^) N x (H^~ 1 ) NxN ) and 
a subsequence (which we still denote by {(p n , u", E n )} n ^), such that, for all p € N, we 
have 

(x P P n , X P u n , x P E n ) -> (x P P, X P u, x P ^) (5.13) 

in C([0,p];ijf " 1 x (H^) 1 * x (iff - 1 ) 7Vx7V ). In particular, this implies that (p n ,u n ,S n ) 
tends to (p,u,E) in D'(M + x R^) as n — >• 00. Furthermore, according to the uniform 
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bounds in the previous subsection, we deduce that (p, u, E) belongs to 

HL 1 (r+;B% +1 >% x (i?f +1 ) w x (B% +1 'f) N * N ) , 

and belongs to C^(R + ; B^' 1 ) x (C* (R+; 5^))^ x (C^ (R + ; B^~ 1 )) NxN (and also be- 
longs to CI(]R + ; Sf - 1 x (Bf- 2 ) N x (flf -y xAr ) if JV > 3). And, obviously, we have the 

bounds provided by (^p n ) for this solution. More important, we can claim that the whole 

— i 
sequence is a Cauchy sequence in 33 2 as the initial data is such small that 2CT7 < ^. 



Indeed, if we denote 

U n = \\{p n -p n - 1 ,u n -u n - 1 ,E n -E n - 1 ) 



N i 
25 ~2~ 



then applying Proposition 12,21 and following the estimates in the previous subsection, we 
can deduce that U n satisfies 

U n (t) < a n +-U n ~ 1 (t), for all t G R+, 
where a n > 0, obtained from the initial data, satisfies ^ ngFJ a n < 00. This inequality 

AT 

easily ensure that the sequence {(p n , u n , E n )} n& fq is a Cauchy sequence in 

Next, we need to prove that (p,u,E) obtained above solves (jl.ip . To this end, we first 
observe that, according to ([37 



(5.14) 



'd t p n+1 + u n • Vp n+1 + divu n+1 = -p n divu n , 
d t u n+l + u n • Vu n+1 - „4u n+1 + Vp n+1 - div£ n+1 

= K(p n )Vp n + j^Au n + div((E n )(E n ) T ), 
d t E n+1 + u n ■ VE n+1 + Vu" = Vu n E n+1 . 

Hence, the only problem now is to pass to the limit in D'(R + x R^) for each term in 
(|5.14p . especially for those nonlinear terms. Let 9 £ C^°(IR + x R") and p £ N be such 
that supp 9 C [0,p] x B(0,p). For the convergence of ^ n Au n , we write 



1 + p n l + p J l + p r ' 

*pP n MJ] Au 

1 + Xp p n 1 + XppJ p 



+ 



n ■ N 

Since 9 1 ^ w is uniformly bounded in L°°(B~) C L°°(H~) and Xp u ™ tends to Xp u m 

■ N — 3 ■ iV— 3 

C([0,p];i? 2 ) ; the first term in the above identity tends to in C([0,p];H 2 ) 5 while 
1 +t P r ," tends to -t+ff _ in C([0,p\; H^" 1 ) by Lemma [2.1( which implies that the third 

■ N — 3 

term also tends to in C([0,p]; H~2~). The convergence of other nonlinear terms can be 
treated similarly. 

Finally, we now prove the continuity of p and E in time. 

~ jV -1 N 

Lemma 5.2. Let (p, u, E) be a solution to (jl.ip . Then p and E are continuous in B~~ '~ , 
andu belongs to C(R + ; B^" 1 ). 
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Proof. To prove this, we follow the argument in [3J. Indeed, the continuity of u is straight- 
forward, because u satisfies 



0, 



u = -u • Vu + Au - Vp - (t^) Au ~ K (P) V P + E jk V Xj E ik , 



and the right-hand side belongs to L l (B^ *) + L 2 {B^ l ) due to the facts that p G 
L 2 (ST - 1 ) and E G L 2 (5t 

To prove p,E £ CQR+jBt- 1 ), we notice that, p ,£ G Bt" 1 , p.fie flf" 1 ) 
and dtp, dtE G L 2 (M + ; Thus, it remains to prove the continuity in time of p, E in 

B~ . To this end, we apply the operator A q to the first equation of (jl.ip to yield 

dtAgp = -A,(u • Vp) - AA q d - A 3 (pdivu). (5.15) 

Obviously, for fixed q, the right-hand side belongs to L 1 (M + ;L 2 ). Hence, each A q p is 
continuous in time with values in L (thus in B~). Now, applying an energy method to 
(|5.15p . thanks to Lemma 12.21 we obtain 

\j t \\^iP\\l 2 ^ Cll NpIIl 2 (V" 9 ^ \\p\\ b k l|u|l B f +i + l|AA 9 d|| L 2 + ||A g (pdivu)||^ 
Integrating the above inequality with respect to time on the interval [ti,^]? we obtain 

2«T||A 9 p(t 2 )|| L2 <2^\\A qP (t 1 )\\ L ,+C ^ {a q {r)\\p{r)\\ B n ||u(r)|| fl ^ +1 
+ 2^f+i)||A 3 d(r)|| L2 +2«T||A 3 (pdivu)(T)|| i2 )dT 

N -> N i-i -1 AT 

Since p G L°°(i?~2~), u G Z/(i?~2~ +1 ), and pdivu G L (2?~2~), we eventually obtain 
l|p(*2)H B f < \\p(ti)\\ B N + (1 + ||p|l LOO(B N } ) J t 2 ||u(r)|| sf +1 dr + jf 2 ||pdivu(r)|| B ^dr, 

which implies that p belongs to C(M + ;S~). 

Similarly, we can prove that E also belongs to C(M + ; B~). 

This finishes our proof. □ 



6. Appendix 

For the completeness of the presentation, we give in this appendix the proof of two 
fundamental lemmas concerning the divergence and curl of the deformation gradient in 
the viscoelasticity system (jl.ip . 

The first one we are going to state now is the following lemma (cf. Proposition 3.1 in 

US]). 

Lemma 6.1. Assume that div(poF,}) = and (p, u, F) is the solution of the system (jl.ip . 
Then the following identity 

div(pF T ) = (6.1) 

holds for all time t > 0. 
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Proof. First, we transpose the third equation in (jl.lj) and apply the divergence operator 
to the resulting equation to yield 



<),{<), F,i) ■ u • V(a,. ( F J ,) 



dxkdxj 



Fjfci- 



(6.2) 



Multiply the first equation in (jl.lj) by d Xj Fji, multiply (j6.2p by p, and summing them 
together, we obtain 



jo 



P 



dxkdxj 



F/,-, - pd Xk u k d Xj Fji. 



(6.3) 



On the other hand, we differentiate the first equation in (jl.ip with respect to Xj to yield 



dt(d x ,p) + u • V{d x .p) + d X] u k d Xk p + 



<9 2 u fc 



p + d Xj pd Xk u k = 0. 



•4) 



dx k dxj 

Multiplying (j6.4p by Fjj, multiplying the third equation in (jl.lj) by c^p, and summing 
them together, we obtain 

d 2 u k 



d t (d Xj pFji) + u • V(d Xj pFji 



-P 



dx k dxj 



■ 31 



d Xj pd Xk u k F 



3i- 



Adding (|6.3p and (|6.5p together yields 

<9 t (div(pF T )) + div(u ® div(pF T )) = 0. 
If (p, u, F) is sufficiently smooth, we multiply (I6.6P by div(/9F T ), we get 



(6.5) 
(6.6) 



div(pF" 



+ div ( u 



div(pF" 



1 



divu 



div(pF _ 



Integrating the above identity with respect to x over R , we obtain 



d 



div(pF n 



2 

L 1 



divu 



div(pF _ 



dx 



< — llVullroo 

- 2 1| Hi 



div(pF n 



2 

L 2 



which, by Gronwall's inequality, implies that, for all t > 



dw(pF T )(t) 



< 



div(/) F c 



L-' 



1 ft 



f*\\Vu(T)\\ L ~dT 



Hence, if <Hv(poFq) = 0, the above inequality will gives lldiv^F 7 )^ = for all t > 0, 
which implies that div(pF T ) = for all t > 0. 

This finishes the proof. □ 

Another hidden, but important, property of the viscoelastic fluids system (jl.lj) is con- 
cerned with the curl of the deformation gradient (for the incompressible case, see |13jll6j). 
Actually, the following lemma says that the curl of the deformation gradient is of higher 
order. 

Lemma 6.2. Assume that (jl.lcp is satisfied and (u, F) is the solution of the system (jl.ip . 
Then the following identity 

F lk ViFij = FijVtfiu (6.7) 
holds for all time t > if it initially satisfies (|6.7p . 
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Proof. First, we establish the evolution equation for the equality F^V/Fy — F/jV/Fj^. In- 
deed, by the equation ([l.lcp . we can get 

dtViFij + u • VVjFj,- + Vzu • VFjj = V m u;V/F mj + V;V m UjF„y . 

Thus, 

Flk(dtViFij + u • VVjFjj) + F/fcV/u • VF^ = F^.V m UjV;F mj ' + F Zfc V z V mUi F mi . (6.8) 
Also, from ([l.lcp . we obtain 

VjFjj^Fifc + u • VF ;fc ) = ViF i3 -V m UjF mJfc . (6.9) 
Now, adding ([6.8[) and (|6.9[) . we deduce that 

d t (F lk ViFij) + u • V(Fj[fcVjF»j) = -Fi k Vi\i ■ VF {j + F ifc V m UiV z F mi 

+ F ;fc V ; V m UiF mj + V«F ij V m u z F mfc (6.10) 
= F/fcV m UjVzF m j + F; fc V;V m UjF m j. 
Here, we used the identity which is derived by interchanging the roles of indices I and m: 
FzfcVzU • VF y - = F; fc VzU m V m Fjj = ViFijV m UiF mk . 

Similarly, one has 

dt(FijS/iF ik ) + u • V(F{jVjFifc) = FijV m UiViF mk + FyVjV m UjF mfc . (6.11) 

Subtracting (|6.1ip from (|6.10|) yields 

c^(Fi fc VjFy - FyVjF ifc ) + u • V(F /fc V z F ij - - F^V/F^) 

= V m Uj(F; fc ViF m j - F;jV;F mfc ) + VjV m Ui(F m jFj fc - F mfe Fy). 

Due to the fact 

V/V m Uj = V m V z Ui 

in the sense of distributions, we have, again by interchanging the roles of indices I and m, 

V;V m Uj(F m jF; fc — F mk F[j) = V;V m UjF mJ F; fc — V;V m UjF. mfc F;j 

= V;V m UjF m jF/fe - \7 m ViUiFi k F m j 
= (V/V m Uj - V m ViUi)Fi k F mj = 0. 
From this identity, equation (|6.12j) can be simplified as 

dt{Fi k ViFij - FyVjFifc) + u • V(F lk ViFij - F Zj V/F ifc ) ^ 
= V m Ui(Fi k ViF m j - FijViF mk ). 
Multiplying dSHSD by FjfcVjFy- - FyVjFifc, we get 

dt\Fi k ViFij - F[jViF ik \ 2 + u • X7\Fi k ViFij - F[jViF ik \ 2 

= 2(F Zfe VzFy - Fy VzF ifc )V m Uj(F/ fe VzF m j - F^VjF mfe ) (6.14) 
< 2||Vu|| Loo(R3) ^ 2 , 
where Af is defined as 

Ai = maxflFjfcVjFy - F y V z F ifc | 2 }. 

Hence, (|6.14j) implies 

dtAl + u- VM < 2||Vu|| L cx>( M 3)A4. (6.15) 
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On the other hand, the characteristics of dtf + u • V/ = is given by 

^- s X(s) = u(s,X(s)), X(t) = x. 
Hence, (|6.14p can be rewritten as 

^ < B(t, y)U, U(0, y) = M (y), (6.16) 

where 

U(t,y)=M(t,X(t,x)), B(t,y) = 2\\Vu\\ Lxim (t,X(t,y)). 
The differential inequality (|6.16p implies that 



Hence, 



U(t,y) < U(0)exp(J B(s,y)ds 
M(t,x) < .M(0)exp ^ 2||Vu|| L oo (IR 3)(s)(is 



Hence, if Ai(0) = 0, then M(t) = for all t > 0, and the proof of the lemma is complete. 

□ 



Using F = I + E, (|6T7|) means 

v k E i:j + /•:,/, v, /.•;, = VjE ik + /-. v, /•:,/,. (6.17) 

According to (|6.17p , it is natural to assume that the initial condition of E in the viscoelastic 
fluids system (|1.3j) should satisfy the compatibility condition 

V fe £((% + E(0)i k ViE(0)ij = VjE(0) lk + £?(0) {j V,£?(0) ifc . (6.18) 
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